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Purification of quantum trajectories 

Hans Maassen Burkhard Kummerer 

Abstract: We prove that the quantum trajectory of repeated perfect mea- 
surement on a finite quantum system either asymptotically purifies, or hits 
upon a family of 'dark' subspaces, where the time evolution is unitary. 



1. Introduction 



A key concept in the modern theory of open quantum systems is the notion of 
indirect measurement as introduced by Kraus [Kra] . An indirect measurement on 
a quantum system is a (direct) measurement of some quantity in its environment, 
made after some interaction with the system has taken place. 

When we make such a measurement, our description of the quantum system changes 
in two ways: we account for the flow of time by a unitary transformation (following 
Schrodinger), and we update our knowledge of the system by conditioning on the 
measurement outcome (following von Neumann). If we then repeat the indirect 
measurement indefinitely, we obtain a chain of random outcomes. In the course of 
time we may keep record of the updated density matrix O t , which at time t reflects 
our best estimate of all observable quantities of the quantum system, given the 
observations made up to that time. This information can in its turn be used to 
| predict later measurements outcomes. The stochastic process 9 t of updated states, 

is the quantum trajectory associated to the repeated measurement process. 
By taking the limit of continuous time, we arrive at the modern models of con- 
tinuous observation: quantum trajectories in continuous time satisfying stochastic 
Schrodinger equations [Dav], [Gis], [Car], [BGM]. These models are employed with 
great success for calculations and computer simulations of laboratory experiments 
C3 such as photon counting and homodyne field detection. 

In this paper we consider the question, what happens to the quantum trajectory 
at large times. We do so only for the case of discrete time, not a serious restriction 
indeed, since asymptotic behaviour remains basically unaltered in the continuous 
time limit. 

We focus on the case of perfect measurement, i.e. the situation where no information 
flows into the system, and all information which leaks out is indeed observed. In 
classical probability such repeated perfect measurement would lead to a further and 
further narrowing of the distribution of the system, until it either becomes pure, i.e. 
an atomic measure, or it remains spread out over some area, thus leaving a certain 
amount of information 'in the dark' forever. Using a fundamental inequality of 
Nielsen [Nie] we prove that in quantum mechanics the situation is quite comparable: 
the density matrix tends to purify, until it hits upon some family of 'dark' subspaces, 
if such exist, i.e. spaces from which no information can leak out. A crucial difference 
with the classical case is, however, that even after all available information has been 
extracted by observation, the state continues to move about in a random fashion 
between the 'dark' subspaces, thus continuing to produce 'quantum noise'. 
The structure of this paper is as follows. In Section 2 we introduce quantum mea- 
surement on a finite system, in particular Kraus measurement. In Section 3 repeated 
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measurement and the quantum trajectory are introduced, and in Section 4 we prove 
our main result. Some typical examples of dark subspaces are given in Section 5. 

2. A single measurement 

Let A be the algebra of all complex dxd matrices. By S we denote the space of d x d 
density matrices, i.e. positive matrices of trace 1. We think of A as the observable 
algebra of some finite quantum system, and of S as the associated state space. 
A measurement on this quantum system is an operation which results in the ex- 
traction of information while possibly changing its state. Before the measurement 
the system is described by a prior state 9 £ S, and afterwards we obtain a piece 
of information, say an outcome i & {1,2, ... , k}, and the system reaches some new 
(or posterior) state 6'f 

(9 — ► (i, «0 . 

Now, a probabilistic theory, rather than predicting the outcome i, gives a probability 
distribution (m, TT2, ■ ■ ■ , TTfe) on the possible outcomes. Let 

Ti : 6 >-> 7Tj#- , (i = l,...,k). (1) 

Then the operations Tj, which must be completely positive, code for the prob- 
abilities ni = tr(Tj#) of the possible outcomes, as well as for the posterior states 
9\ = Ti9 /tr(Tj0), conditioned on these outcomes. The fc-tuple (Ti, . . . , Tfc) describes 
the quantum measurement completely. Its mean effect on the system, averaged over 
all possible outcomes, is given by the trace-preserving map 

k k 
T : 9 ^^TTiO'i = ^T t 9 . 

i=l i=l 

Example 1: von Neumann measurement. 

Let pi,P2, ■ ■ -Pk be mutually orthogonal projections in A adding up to 1, and let 
a G A be a self- adjoint matrix whose cigcnspaces are the ranges of the pi. Then 
according to von Neumann's projection postulate a measurement of a is obtained 
by choosing for Tj the operation 

T,{9)=p,9p, . 
Example 2: Kraus measurement. 

The following indirect measurement procedure was introduced by Karl Kraus [Kra] . 
It is contains von Neumann's measurement as an ingredient, but is considerably 
more flexible and realistic. 

Our quantum system A in the state 9 is brought into contact with a second system, 
called the 'ancilla', which is described by a matrix algebra B in the state 0. The 
two systems interact for a while under Schrodinger's evolution, which results in 
a rotation over a unitary u € B ® A. Then the ancilla is decoupled again, and 
is subjected to a von Neumann measurement given by the orthogonal projections 
Pi,...,Pk € B. The outcome of this measurent contains information about the 
system, since system and ancilla have become correlated during their interaction. 
In order to assess this information, let us consider an event in our quantum system, 
described by a projection q e A. Since each of the projections Pi ® 1 commutes with 
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1 ® q, the events of seeing outcome i and then the occurrence of q are compatible, 
so according to von Neumann we may express the probability for both of them to 
happen as: 

P[outcome i and then event q] = tr <g> tr^(u(/3 ® {Pi <8> ?)) • 

Therefore the following conditional probability makes physical sense. 

tr®tr ((«(/? ®0)u*)(pi<8i«)) 



P [event g| outcome i] 



tr®tr((u(/3®0)u*)(pj ® 1)) 



This expression, which describes the posterior probability of any event q £ A, can 
be considered as the posterior state of our quantum system, conditioned on the 
measurement of an outcome i on the ancilla, even when no event q is subsequently 
measured. As above, let us therefore call this state Q\. We then have 



tr(T, 

where T;B takes the form 



Tj0 = te®id((u(/?®0)u*)(pi<8il)) . 



Here, id denotes the identity map S — > 5. 

The expression for takes a simple form in the case which will interest us here, 
namely when the following three conditions are satisfied: 

(i) B consists of all k x k- matrices for some k; 

(ii) the orthogonal projections Pi G B are one-dimensional (say Pi is the matrix 
with i-th diagonal entry 1, and all other entries 0); 

(iii) (5 is a pure state (say with state vector (/3i, . . . , /3k) € C fe ). 

These conditions have the following physical interpretations. 

(i) The ancilla is purely quantummechanical; 

(ii) the measurement discriminates maximally; 

(iii) no new information is fed into the system. 

If these conditions are satisfied, u can be written as a k x k matrix (uij) of d x d 
matrices, and may be written 

TiB = ai 9a* , (2) 

where 

k 

a * = H fo^i ■ 

We note that, by construction, 

i— 1 i— 1 j — 1 j' — 1 

This basic rule expresses the preservation of the trace by T. 
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Definition 1 By a perfect measurement on A we shall mean a k-tuple (Ti, . . . , Tk) 
of operations on S, where Tfi is of the form QiOa* with X^j=i a \ a i = 1- 

Mathematically speaking, the measurement (Ti, . . . , Tk) is perfect iff the Stincspring 
decomposition of each Tj consists of a single term. 

We note that every perfect Kraus measurement is a perfect measurement in the 
above sense, and that every perfect measurement can be obtained as the result of 
a perfect Kraus measurement. 

3. Repeated measurement 

By repeating a measurement on the quantum system A indefinitely, we obtain a 
Markov chain with values in the state space S. This is the quantum trajectory 
which we study in this paper. 

Let fi be the space of infinite outcome sequences ui = W2, W3, . . .}, with uij € 
{1, . . . , k}, and let for m e N and ii, ■ ■ ■ ,i m & {1, . . . , k} the cylinder set A il ^,^ im C 
n be given by 

K lt ...,i m ■= W e fi| uj-i = i u . .. ,u> m = i m } . 

Denote by S m the Boolean algebra generated by these cylinder sets, and by S the 
a- algebra generated by all these S m . Let Ti, . . . , Tk be as in Section 2. 
Then for every initial state 80 on A there exists a unique probability measure Fg 
on (ft, E) satisfying 

Poo(Ai 1 ,...,i m )=tr(T ira o...oT <1 (e )) • 

Indeed, according to the Kolmogorov-Daniell reconstruction theorem we only need 
to check consistency: since T = ^ i=1 Tj preserves the trace, 

k k 

P <>« ( A h,.,im,i) = H tr (^ °T tm o---oT ll (0 O )) = tr(T o T im o • • • o T h (9 )) 

i=l i=l 

= tr(T im o...oT il (fl )) =Po (Ai 1 ,...,i m ) • 



On the probability space (fl, S, Pg ) we now define the quantum trajectory (@n) ne f$ 
as the sequence of random variables given by 

B„ : ll — > o : uj 



tr(T Um o.-oT ul (fl )) ' 

We note that 9„ is S„-measurable. The density matrix 6 n (w) describes the state 
of the system at time n under the condition that the outcomes w\,...,w n have been 
seen. 

The quantum trajectory (6 n )„ e p^ is a Markov chain with transitions 

d ^ d ' t = t^rT) with P robabilit y ^(T t e). (3) 
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4. Purification 

In a perfect measurement, when T, is of the form 8 i— > ai6a*, a pure prior state 
9 — IV'XV'I leads to a pure posterior state: 

, ai\ip)(ip\a* diip 
d i = -T, — * — 7T Hwwl , where Vi - 



(ip, a*a t ip) ' \\aiip\\ 

Hence in the above Markov chain the pure states form a closed set. Experience with 
quantum trajectories leads one to believe that in many cases even more is true: along 
a typical trajectory the density matrix tends to purify: its spectrum approaches the 
set {0, 1}. In Markov chain jargon: the pure states form an asymptotically stable 
set. 

There is, however, an obvious counterexample to this statement in general. If every 
a, is proportional to a unitary, say a; = ^/XiUi with u*Ui = 1, then 

UiVU; <~ 9 , 



1 tv( ai 6a*) 

where ~ denotes unitary equivalence. So in this case the eigenvalues of the density 
matrix remain unchanged along the trajectory: pure states remain pure and mixed 
states remain mixed with unchanging weights. In this section we shall show that 
in dimension 2 this is actually the only exception. (Cf. Corollary 2.) In higher 
dimensions the situation is more complicated: if the state does not purify, the ai 
must be proportional to unitaries on a certain collection of 'dark' subspaces, which 
they must map into each other. (Cf. Corollary 8.) 

In order to study purification we shall consider the moments of O n . By the m-th 
moment of a density matrix 9 € S we mean tr (0 m ). We note that two states 9 and 
p are unitarily equivalent iff all their moments are equal. In dimension d equality 
of the moments m = 1, . . . , d suffices. 

Definition 2 We say that the quantum trajectory (6 n (w)) ne ^ purifies when 
W meN-- V* tr(enH m )=l. 

Note that the only density matrices p satisfying tr(p m ) — 1 are one-dimensional 
projections, the density matrices of pure states. In fact, it suffices that the second 
moment be equal to 1. 

We now state our main result concerning repeated perfect measurement. 

Theorem 1 Let (©n) ng p^ be the Markov chain with initial state 9 and transition 
probabilities (3). Then one of the following alternatives holds. 

(i) The paths of (Q n ) ne ^q (the quantum trajectories) purify with probability 1, or: 
(ii) there exists a projection p e A of dimension at least two such that 

Vie{i,...,fe}3A i >o : pa*(HP = ^iP ■ 

Condition (ii) says that a, is proportional to an isometry in restriction to the range 
of p. Note that this condition trivially holds if p is one-dimensional. 

Corollary 2 In dimension d — 2 the quantum trajectory of a repeated perfect mea- 
surement either purifies with probability 1, or all the ai 's are proportional to uni- 
taries. 
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If the a, are all proportional to unitaries, the coupling to the environment is essen- 
tially commutative in the sense of [KiiM]. 

Our proof starts from an inequality of Michael Nielsen [Nic] to the effect that for 
all m e N and all states 9: 

k 

5>* tr ($r) >tr(6n, 

i=l 

where 

7Ti := tr(ai#a*) and ^ := 



tr(aj#a|) 

Nielsen's inequality says that the expected m-th moment of the posterior state is 
as least as large as the m-th moment of the prior state. In terms of the associated 
Markov chain we may express this inequality as 

V m ,„ e N : E(tr(6™ +1 )|s n ) >tr(9") , 

i.e. the moments M n m ^ := tr (@™)„ e j^ are submartingales. Clearly all moments 
take values in [0, 1]. Therefore, by the martingale convergence theorem they must 
converge almost surely to some random variables M^ m \ 

This suggests the following line of proof for our theorem: Since the moments con- 
verge, the eigenvalues of (9n)„ e p^ must converge. Hence along a single trajectory 
the states eventually become unitarily equivalent, i.e. eventually 

w n f a <k®n(u)a* 
Vj : 0„(w) - 



tr(oi8„(w)a*) ' 

But this seems to imply that either 0„ purifies almost surely, or the a^'s are unitary 
on the support of 6„. 

In the following proof of Theorem 1 we shall make this suggestion mathematically 
precise. 

Lemma 3 In the situation of Theorem 1 one of the following alternatives holds, 
(i) For all m <G N: lim tr (9™) = 1 almost surely; 

n— >oo 

(ii) there exists a mixed state p e S such that 

Vi=i,...,fe3 Ai >o : aipa* <~ \ t p . 



k / / / n * \ 'tn\ \ 2 



tv(a,i0a*) 



tr{6 m ) 



Proof. For each m € N we consider the continuous function 

k . 
6 m : S -» [0, oo) : 9 ^ ^ tr(a,(9a*) ( tr 
i=i ^ 

Then, using (2) and (3), 

Since (M„ ) n£ j^ is a positive submartingale bounded by 1, its increments must 
be square summable: 

oo 

v meN : EE(Me n )) <i. 

n=0 
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In particular 



lim VE(<5 m (6„)) = 0. (4) 



n^oo 

m— 1 



Now let us assume that (i) is not the case, i.e. for some (and hence for all) m > 2 
the expectation E(M^) — : /z m is strictly less than 1. For any jigN consider the 
event 

A n :={c^|m( 2 ><^±±} . 
Then, since is increasing in n, we have for all n € N: 

E(M^) <E(M< 2 >) =/i 2 < 1 . 

Therefore for all n € N, 



^2 > E(M( 2 )-l [M ( 2)>ii ^ i] ) 



> 



M2 + 1 

2 

^2 + 1 



Mi 2 ) 



> 



M2 + 1 



(i-P(A,)) 



so that 



F(A„) > - 



1 - ^2 
M2 



(5) 



On the other hand, A„ is £„-measureable and therefore it is a union of sets of the 
from Aj lv .. ; j n . Since 6„ is S„-measurcable, 9„ is constant on such sets; let us call 
the constant Q n (h, ■ ■ ■ , in)- We have the following inequality: 



1 ( d \ 1 d 

n ' A i;1 inC A„ \m=l / F ^"' ) m=l 



e„)) . 



On the left hand side we have an average of numbers which are each of the 
form J2m=i $m(®n(ii, ■ ■ ■ , in)), hence we can choose . . . ,i n ) such that p n := 
Q n (ii, . . . , i n ) satisfies, by (5), 

M2 + 1 



E MpJ < ———— V E(j ro (e n )) . 

— i A*2 — t — . 

m— 1 m— 1 

Since Aj^...^ C A n , the sequence (Pn)„ e N nes en thely in the compact set 

{ees tr(e 2 )<^±l 

Let p be a cluster point of this sequence. Then, since E(<$ m (0„)) tends to as 
n — > oo, and S m is continuous, we may conclude that for m = 1, . . . , d: 



5 m (p) = 0, and tr(p 2 ) < ^±1 < 1 . 
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So p is a mixed state, and, by the definition of S m , 

"(»«»'••> ( tr ((«)") - tt( "™ , ) 2 = 

for all m = 1, 2, 3, . . . , d and all z = 1, . . . , k. Therefore either tr(aipa*) = 0, i.e. 
aipa* = 0, proving our statement (ii) with A,; = 0; or tr(a,ipa*) > 0, in which case 
p\ := dipa* /tv{aipa*) and p itself have the same moments of orders m = 1, 2, . . . , d, 
so that they are unitarily equivalent. This proves (ii). 

□ 

From Lemma 3 to Theorem 1 is an exercise in linear algebra: 

Lemma 4 Let ai,...,a,k € Md be such that Y^l=i a * a i = 1 Suppose that there 
exists a density matrix p € Md such that for i = 1, . . . , k 

aipa* ~ Xip . 

Let p denote the support of p. Then for all i = 1, . . . , k: 

pa*aip = Xip . 



Proof. Let us define, for a nonnegative matrix x, the positiue determinant det pos (a;) 
to be the product of all its strictly positive eigenvalues (counted with their multi- 
plicities). Then, if p denotes the support projection of x, we have the implication 

detp OS (x) = dctpos(Ap) => tr(xp) > tr(Ap) (6) 

with equality iff x = Xp. (This follows from the fact that the sum of a set of positive 
numbers with given product is minimal iff these numbers are equal.) 
Now let p be the support of p as in the Lemma. Let Vi^/pa*aip denote the polar 
decomposition of aip. Then we have by assumption, 

det pos (Xip) = dct pos (aipa*) 

= det pos (aipppa*) 

= det pos (vi^pa*aipp^pa*aipv*) 

= detp OS ( y/ pa* a { pp^J pa* Oip) 

= det pos (pa*aip)det pos (p) . 

Now, since det pos (Xip) — det pos (Aip) • det pos (/9) and det pos (/5) > 0, it follows that 

det pos (Ai_p) = dct pos (pa*aip) . 

By the implication (6) we may conclude that 

tv(pa*aip) > trX l p) . (7) 

On the other hand, 

fe k k // fc \\ 

Xl = tr ( A ^) = tT ( a *P a *t) = tr I P I Y ata *i J = tvp = 1 ' 
i=i i=i i=i \ \i=i ) ) 
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where in the second equality sign the assumption was used again. Then, by (7), 

k k / k \ 

trp = tr(pa*aip) > tr(Xip) — I Ai ) tip = tip . 

i=l i=l \i=l / 

So apparently, in this chain, we have equality. But then, since equality is reached 
in (6), we find that 

pa*aip = Xip . 



□ 



5. Dark subspaces 



By considering more than one step at a time the following stronger conclusion can 
be drawn. 

Corollary 5 In the situation of Theorem 1, either the quantum trajectory purifies 
with probability 1 or there exists a projection p of dimension at least 2 such that for 
all I £ N and all i\, . . . ,i\ there is X it it > with 

pa* t ■■■a*a n ■ ■ ■ a h p = A, , /» . (8) 

We shall call a projection p satisfying (8) a dark projection, and its range a dark 
subspace. 

Let p be a dark projection, and let Vi^Jpa*ai = \f\~iVip be the polar decomposition 
of aip. Then the projection p\ :— Vipv* satisfies: 



X iPi a n ■ ■ ■ a l m a i m ■ ■ ■ a hPi = KiviPv*)®*! ■ ■ ■ a * m a i m ■ ■ ■ a h (viPV*) 

= Vipa^a*^ ■ ■ ■ a* m a im ■ ■ ■ a^aipv* 

\,ii r ..,! m ' Pi • 

Hence if p is dark, and Ai ^ then also p\ is dark with constants 

\x,...,i m — \,h,...,i m /^i ■ 

We conclude that asymptotically the quantum trajectory performs a random walk 
between dark subspaces of the same dimension, with transition probabilities p — ► 
p\ equal to Ai, the scalar value in pa*aip = \p. In the trivial case that the dimension 
of p is 1, purification has occurred. 

Inspection of the ai should reveal the existence of nontrivial dark subspaces. If none 
exist, then purification is certain. 

We end this Section with two examples where nontrivial dark subspaces occur. 

Example 1. Let d = l-e and let Hi, ... ,Hi be mutually orthogonal e-dimensional 
subspaces of H = C d . Let (iTij) be an I x I matrix of transition probabilities. Define 
fly e A by 

On . — \/'KiiVii , 



where the maps Vij : Hi — > Hj arc isometric. Then the matrices aij, i,j = 1, . . . , 
define a perfect measurement whose dark subspaces are Hi , . . . , Hi . 
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Example 2. 

The following example makes clear that nontrivial dark subspaces need not be 
orthogonal. 

Let Ti — C 2 <£> V, where V is some finite dimensional Hilbert space, and for i = 
1, . . . , k let di := bi ® u,, where the 2 x 2-matrices bi satisfy the usual equality 

k 

and the U{ are unitaries I? — > T>. Suppose that the bi are not all proportional to 
unitaries. Then the quantum trajectory defined by the otj has dark subspaces ^®T>, 
with ip running through the unit vectors in C 2 . Physically this example describes a 
pair of systems without any interaction between them, one of which is coupled to 
the environment in an essentially commutative way, whereas the other purifies. 



References 

[BGM] L. Boutcn, M. Gut,a, H. Maassen: Stochastics Schrddinger equations, 
Optics and Spectroscopy 94 (2003) 911-919. 

[Car] H.J. Carmichael: An Open systems Approach to Quantum Optics, 

Springer- Verlag, Berlin Heidelberg, New York 1993. 

[Dav] E.B. Davies: Quantum Theory of Open Systems, Academic Press, Lon- 

don, New York, San Francisco 1976. 

[Gis] N. Gisin: Quantum measurements and stochastic processes, Phys. Rev. 

Lett. 52 (1984) 1657-1660. 

[Kra] Karl Kraus: States, effects and operations: fundamental notions of 

quantum theory, Lecture Notes in Physics, Vol. 190. Springer- Verlag, 
Berlin 1983. 

[KiiM] Burkhard Kummerer, Hans Maassen: The essentially commutative di- 

lations of dynamical semigroups on M n , Commun. Math. Physics 109 
(1987) 1-27. 

[Nie] Michael Nielsen: Characterizing mixing and measurement in quantum 

mechanics, Physical Review A 63 (2000) art. nr. 022114. 

Radboud University Nijmegen Technische Universitdt Darmstadt 

e-mail: maassenOmath.ru.nl e-mail: kuemmererOmathematik.tu-darmstadt.de 



imsart-lnms ver. 2005/02/28 file: Maassen.tex date: February 1, 2008 



